Half-Heusler Compounds as a New Class of Three-Dimensional Topological Insulators 
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Using first-principles calculations within density functional theory, we explore the feasibility of 
converting ternary half-Heusler compounds into a new class of three-dimensional topological insu- 
lators (3DTI). We demonstrate that the electronic structure of unstrained LaPtBi as a prototype 
system exhibits distinct band-inversion feature. The 3DTI phase is realized by applying a uni- 
axial strain along the [001] direction, which opens a bandgap while preserving the inverted band 
order. A definitive proof of the strained LaPtBi as a 3DTI is provided by directly calculating the 
topological Z2 invariants in systems without inversion symmetry. We discuss the implications of 
the present study to other half-Heusler compounds as 3DTI, which, together with the magnetic 
and superconducting properties of these materials, may provide a rich platform for novel quantum 
phenomena. 
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Recent years have seen a surge of interest in a new class 
of materials called topological insulators [THSJ- These 
materials are distinguished from ordinary insulators by 
nontrivial topological invariants associated with the bulk 
electronic structure [IH5] . The existence of the topologi- 
cal invariants dictates that the excitation gap must van- 
ish at the boundaries of a topological insulator, result- 
ing in the formation of robust metallic surface states. A 
number of spectacular quantum phenomena have been 
predicted when the surface states are under the influence 
of magnetism and superconductivity [Tllllj . To fully ex- 
plore these phenomena thus demands great versatility 
from the host material. However, so far the experimen- 
tal realizations of topological insulators arc limited to a 
few classes of simple materials, including HgTe quantum 
well [12 H3], Bii-zSbz alloy [TS], and tetradymite 
semiconductors such as Bi 2 Se3, Bi 2 Te3, and Sb2Te3 [TBl - 
[18] . Realizing the necessary conditions for the predicted 
phenomena in these materials can be difficult. 

The search for topological insulators has greatly bene- 
fited from the topological band theory. It has been shown 
that for all known topological insulators, in addition to 
the strong spin-orbit coupling, their electronic structure 
can be characterized by a band-inversion which involves 
the switching of bands with opposite parity around the 
Fermi level [T5J HH US] • This is very similar to the quan- 
tum Hall effect in which two bands are allowed to ex- 
change their Chern numbers only when they come into 
contact with each other [TH] . The above observation sug- 
gests small bandgap semiconductors or semimetals with 
heavy elements as promising candidates as these materi- 
als are likely to develop an inverted band order. 

In this Letter we predict a new class of topological in- 



sulators realized in small bandgap ternary half-Heusler 
compounds [20j . In particular, using first-principles cal- 
culations we demonstrate that LaPtBi as a prototype 
system becomes a strong topological insulator upon the 
application of a uniaxial strain along the [001] direction. 
This result is first discussed using the aforementioned 
band-inversion mechanism, then verified by direct calcu- 
lation of the Z 2 topological invariants from the bulk band 
structure [22 [55]. Remarkably, ternary half-Heusler com- 
pounds already boast an impressive list of the much de- 
sired properties such as magnetism [23] and supercon- 
ductivity [24] • Together with the predicted topological 
properties, these materials provide an exciting platform 
for novel quantum phenomena. 

Ternary half-Heusler compounds have the chemical for- 
mula XYZ, where X and Y are transition or rare earth 
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FIG. 1: (color online) Crystal structure of half-Heusler com- 
pound LaPtBi in the F43m space group. Green spheres at 
(0.5,0.5,0.5) are atom La, dark blue spheres at (0.25,0.25,0.25) 
are atom Pt, and pink spheres at (0,0,0) are atom Bi. 
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FIG. 2: Band structure of LaPtBi with experimental lattice 
constant a=6.83 A. The band order at the T point is Tg, rV, 
and T§ in descendent order of energy. There is a small gap 
between the F7 and F6 states. 



FIG. 3: Band structure of LaPtBi under uniaxial strain with 
constant volume along [001]-direction, reduce the c/a ratio by 
5% in (a), and increase the c/a ratio by 5% in (b). 



metals and Z a heavy element. Figure [T] shows the crys- 
tal structure of LaPtBi, which consists of three interpene- 
trating, face-centered-cubic lattices with Pt sitting at the 
unique site. It can be regarded as a hybrid compound of 
LaBi with the rock-salt structure, and LaPt and PtBi 
with the zinc-blende structure. Unlike tetradymite semi- 
conductors, the spatial inversion symmetry is broken in 
half-Heusler structure. 

When the total valence electron count in a primitive 
unit cell is 18, the half-Heusler compounds are expected 
to have a band gap. However, some compounds display 
a distinctive semimetal behavior with LaPtBi being one 
of the examples [55] • To investigate the band topol- 
ogy, we employ the full-potential linearized augmented 
plane-wave method [2 6) with the local spin density ap- 
proximation for the exchange-correlation potential |27j . 
Fully relativistic band calculations were performed using 
the program package WIEN2K [28]. A converged ground 
state was obtained using 10,000 k points in the first Bril- 
louin zone and K max R,MT = 9, where Rmt represents 
the muffin-tin radius and K mayi the maximum size of the 
reciprocal-lattice vectors. Wave functions and potentials 
inside the atomic sphere are expanded in spherical har- 
monics up to I = 10 and 4, respectively. Spin orbit cou- 
pling are included by a second- variational procedure [55] , 
where states up to 9 Ry above Fermi energy are included 
in the basis expansion, and the relativistic pi/2 correc- 
tions were also considered for 5p, 6p of Pt, and 6p of 
Bi in order to improve the accuracy [29l [30] . The cal- 
culations were performed using the experimental lattice 
constant of 6.83 A [31] . Figure [2] shows the fully rela- 
tivistic energy band structure of LaPtBi. As anticipated 
from experiments [32] . LaPtBi is a semimetal with very 
small electron and hole pockets around the T point. Our 
result is consistent with previous calculations [25 . 

As already pointed out by Oguchi [35], the band 



structure near the Fermi level at the T point is deter- 
mined mainly by PtBi with the zinc-blende structure, 
and the La states participate in the band structure ad- 
ditively. This separation allows us to draw direct com- 
parison with a known topologically nontrivial compound 
HgTe rrjj [II], which is a II- VI material also with the 
zinc-blende structure. Let us focus on the bands at the 
r point close to the Fermi level. Similar to the case of 
HgTe, symmetry analysis shows that the fourfold degen- 
erate r 8 states lies above the twofold degenerate T 7 and 
F 6 states. As discussed by Fu and Kane [TJ, such a band 
inversion is a strong indication that LaPtBi is in a topo- 
logically nontrivial state. To remove the semi-metallic 
behavior, we apply a uniaxial strain along the [001] di- 
rection with constant volume to break the fourfold de- 
generacy of the Tg states. Figure [3] shows the resulting 
band structure. We find that although a local bandgap 
can be opened, the system responds differently to com- 
pression and elongation: Upon compression the material 
remains a semimetal while it becomes an insulator when 
stretched. The inverted band order stays the same. For 
comparison, we also calculated the band structure with 
hydrostatic strain, shown in Fig. [4j A global band gap is 
opened by the hydrostatic strain when the lattice is com- 
pressed while the semi-metallic behavior is retained with 
expansion. However, in the former case the Tq states now 
jumps above the T% states with the Fermi energy lies in 
between. In this situation, the material should be in a 
topologically trivial phase. 

Although the band inversion near the T point is a 
strong indication that LaPtBi under uniaxial strain is 
in a topologically insulating phase, it is not definitive 
because the topological invariant is a global property of 
the entire Brillouin zone. Fu and Kane have proposed a 
parity criterion to identify topological insulators in sys- 
tems with both time-reversal and spatial inversion sym- 
metry [14] . However, it cannot be applied here because of 
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FIG. 4: Band structure of LaPtBi under hydrostatic strain, 
a=a-7%a in (a), a=a+7%a in (b) 



the lack of inversion symmetry in the half-Heusler struc- 
ture. Instead, we calculate the topological invariants di- 
rectly from the bulk band structure. 

We first briefly describe the formalism for a 2D system. 
In the presence of time-reversal symmetry, Kramer's the- 
orem dictates that the energy eigenstates must come 
in pairs. This allows us to enforce the so-called time- 
reversal constraint on the Bloch functions: 



M-*0) = e|«n(fc)) 



(i) 



where |u n (fe)) is the periodic part of the Bloch function, 
and O = e l ^ s yl n K is the time reversal operator with 
S y the spin operator and K the complex conjugation. 
Accordingly, we only need to obtain Bloch functions in 
half of the Brillouin zone, denoted by B + , as those in the 
other half are fixed by Eq. (JlJ. The band topology is 
characterized by the Z 2 invariant, given by |21) 



1 r 

2^ 



dk ■ A(k) 



d 2 kT{k) mod 2, (2) 



where A(k) = i^ n (u n {k)\V ' ku n {k)) is the Berry con- 
nection and JF(k) = x A(k)\ z is the Berry curvature; 
the sum is over occupied bands. A topological insulator 
is characterized by Z 2 = 1 while ordinary insulators have 
Z 2 = 0. The nonzero Z 2 invariant is an obstruction to 
smoothly defining the Bloch functions in B + under the 
time- reversal constraint. 

To numerically perform the integration, we follow the 
recipe by Fukui and Hatsugai 22J. The Bloch func- 
tions \u n (k)) are first obtained on a fe-space mesh in B + . 
The mesh must include the four time-reversal invariant 
fc-points: 0, Gi/2, G 2 /2 and (Gi + G 2 )/2, expressed 
in terms of the reciprocal lattice vectors. After apply- 
ing the time-reversal constraint, next we introduce the 
link variable central to many Berry-phase related calcu- 
lations |33J [M], given by J7 M (fej) = det \\{u n (kj)\u m (kj + 
/x))||, where ft is the unit vector on the mesh, and n and 
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FIG. 5: The n-field configuration for LaPtBi under uniaxial 
strain computed under the time-reversal constraint. The four 
tori are Zq, Zi, Xo and Yo with the shaded area indicating 
half of the area, fi and v are the unit vectors of the fe-space 
mesh. The white and black circles denote n = 1 and — 1, 
respectively, while the blank denotes 0. The Z2 invariant for 
each individual tori is obtained by summing the n-field over 
half of the tori. These read zo = 1, Z\ = 0, xo = 1, and 
yo — 1. The Z2 invariants of the system are 1;(000). 



m run through occupied bands. The finite element ex- 
pressions for A and T are A^ikj) = Imlog U^kj), and 
Ffa) = JmlogU^U^kj + riU-^kj + ^U-^kj), 
where the return value of the complex logarithm function 
is confined to its principal branch (— tt, 7t]. We can then 
insert these expressions into Eq. ([2| to calculate the Z 2 
invariant. Define an integer field n(kj) for each plaque- 
tte: 



n(kj 



^{[A v A^kj 



A/J.A, 



,(k 3 )]~T(k 3 )} , (3) 



where A p is the forward difference operator. The Z 2 
invariant is given by the sum of the n-field in half of the 
Brillouin zone [2"2"] . i.e., Z 2 = 5Z fe . gB + n(kj) mod 2. 

In 3D the topology of the bands are characterized by 
four independent Z 2 invariants @J [S] . These can be com- 
puted by considering six tori in the Brillouin zone. For 
example, the torus Z is spanned by G\ and G 2 with the 
third component fixed at 0, and Zi is obtained by fixing 
the third component at G3/2. The other four tori J^o, 
Xi, Yo, and Yj are defined similarly. For each torus, one 
can calculate the corresponding Z 2 invariant using the 
steps outlined above for 2D systems. Out of the six pos- 
sible Z 2 invariants only four of them are independent. 
Following Ref. [H |S], we use the notation uq; {v\Viv-$) , 
with z/q = (zq + zi) mod 2, v\ = xi, = y\ and 1^3 = Zi, 
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where zq is the Z2 invariant associated with the 2D torus 
Zq. The other Z2 invariants are defined similarly. A 
nonzero vq indicates that the system is a strong topolog- 
ical insulator. 

Figure [5] shows the n-field configuration for LaPtBi 
under uniaxial strain from first-principles calculations. 
The corresponding band structure is shown in Fig. J3jb). 
Note that although the n-field itself depends on the gauge 
choice, their sum over half of the Brillouin zone is gauge- 
invariant module 2. [35J We find that LaPtBi under uni- 
axial strain becomes a strong topological insulator with 
the topological invariants 1;(000). However, under hy- 
drostatic strain [Fig. |4]ja)] , LaPtBi becomes an ordinary 
insulator. Furthermore, although the band structure in 
Fig. [3]Ja) shows a semi-metallic phase, we can still cal- 
culate the Z2 invariants for the bands because a local 
energy gap separates the conduction and valance bands 
throughout the Brillouin zone. In this case, LaPtBi be- 
comes a topological metal with i>o = 1. 

Having firmly established that LaPtBi under uniax- 
ial strain realizes a topological insulating phase, we 
have calculated a number of other half-Heusler ternary 
compounds by first-principles method. We find that 
LuPtSb, ScPtBi, YPdBi, YPtSb have inverted band 
structure, and the 3DTI phase can be induced by uniax- 
ial strain along [001]-dircction. Generally, for many small 
bandgap half-Heusler compounds, the topologically insu- 
lating phase can be realized by a combination of hydro- 
static strain to change the band order and uniaxial strain 
to open an energy gap. Details of the first-principles cal- 
culations will be reported elsewhere [35] . 

In conclusion, we have shown that the 18-electron 
ternary half-Heusler compounds can be tuned into a 
new class of three-dimensional topological insulators via 
proper strain engineering. This is confirmed by first- 
principles calculation of the topological Z 2 invariants in 
systems without inversion symmetry. This quantum na- 
ture, plus other interesting physical properties of these 
materials, such as magnetism [23] and superconductiv- 
ity [24], characterize these materials as an exciting plat- 
form for novel quantum phenomena. 

Note added: After the completion of the bulk of this 
work (see, e.g., the brief announcement in Ref. [2U])> two 
more related studies have appeared |361 137j , confirming 
and expanding the predictions of the present work. 
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